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Abstract A quasi-non-linear fracture mechanics model 
based on beam on elastic foundation theory is applied for 
analysis of the double cantilever beam (DCB) specimen for 
determination of fracture energy of wood. The properties of 
the elastic foundation are chosen so that the perpendicular- 
to-grain tensile strength and fracture energy properties of 
the wood are correctly represented. It is shown that this 
particular choice of foundation stiffness makes a conven¬ 
tional maximum stress failure criterion lead to the same 
solution as the fracture mechanics compliance method. Re¬ 
sults of linear elastic fracture mechanics are obtained as a 
special case by assuming an infinitely large value of the 
perpendicular-to-grain tensile strength. The quasi-non- 
linear fracture mechanics model is compared with other 
models and with results of tests conducted to reveal the 
influence of the geometrical properties of the DCB speci¬ 
men. In addition, the appropriateness of choice of the foun¬ 
dation stiffness is investigated. 

Key words DCB specimen ■ Fracture energy • Beam on 
elastic foundation ■ Foundation stiffness • Quasi-non-linear 
fracture mechanics 


Introduction 

Mode I fracture energy of wood may be determined experi¬ 
mentally using, e.g., small tension specimens 1 or a notched 
specimen in a three-point bending test 2 recording the 
complete load-deflection relationship. A much easier ex¬ 
perimental method, however, is to use a double cantilever 
beam (DCB) specimen as shown in Fig. 1. The DCB speci¬ 
men has been adopted in ASTM D3433 3 for determination 
of the critical energy release rate of bond lines using ad- 
herends of steel, but no major testing standard seems yet to 
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have adopted the DCB specimen for determination of 
the mode I critical energy release rate or the fracture energy 
of wood. 

Using the DCB specimen, however, the critical energy 
release rate is not determined directly by testing, but is 
derived from test results through some theoretical models. 
Thus, the value of the critical energy release rate depends 
on the assumptions adopted in the theoretical model. A 
number of theoretical expressions for determination of the 
critical energy release rate for wooden DCB specimens 
have previously been presented, 1 ’ 4 ' 5 some of which are based 
on beam on elastic foundation (BEF) theory. 

The model presented in this article is likewise based 
on BEF theory. The main new feature in the model pre¬ 
sented here is the definition of the foundation modulus. 
Previous models have associated the foundation modulus 
with the elastic strain in the wood, i.e., a function of the 
perpendicular-to-grain modulus of elasticity (MOE) and 
the dimension of the wood member in the perpendicular-to- 
grain direction, whereas the present model associates the 
foundation modulus with the fracture performance of the 
wood, i.e., perpendicular-to-grain tensile strength and mode 
I fracture energy. Furthermore, the present model is based 
on Timoshenko beam theory, i.e., the shear deformations 
are taken into account, and a finite length of the initially 
uncracked part of the specimen is considered. 


Theory 

A beam on elastic foundation is schematically shown in 
Fig. 2. The deflections and rotations of the beam axis are 
denoted w(x) and 0(x), respectively. Of special interest here 
is w 0 = w(0) and 0 O = 0(0), of which positive directions 
are indicated in the figure. 

For a Timoshenko beam on elastic foundation, 6 w„ and 6 0 
may be written 

P c 3 + qfl _ P q + c 2 a 


( 1 ) 
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Fig. 1. Geometry of double cantilever beam (DCB) specimen 


< 5 P /2 



0 


Fig. 2. Symmetrical half of DCB specimen modeled as beam on elastic 
foundation 


where the definitions of c, c 1; c 2 , and c 3 are given by Eq. 3, 
and P is the applied load as shown in Figs. 1 and 2. 

In addition, the parameters X and rj are defined as 


X = 


Kb 
El ’ 


Kb 

GA 


( 2 ) 


where b is the width, I is moment of inertia, and A s is 
equivalent shear area (I = bh 3 l 12 and A s = 5bh/6 for a 
rectangular cross section), E is the parallel-to-grain modu¬ 
lus of elasticity, G is the shear modulus, and K is the foun¬ 
dation modulus (unit: N/m 3 ). 

The solution to the governing differential equations is 
divided into two cases, and the c constants in Eq. 1 are given 
by 

Case I: X 5® V 477 2 

2 1 [7 , 1 2 1/7 1 

2 4 7 2 4 7 


c — v/2[n 2 sinh 2 v(L — a) — v 2 sin 2 zz(L — a)] 
q = v 2 sin 2 zz(L — a) + zz 2 sinh 2 v(L — a) 
c 2 = 2vz<[acosh v(L — a)sinh v(L — a) 

+ v cos zz(L — a) sin zz(L — a)J 

c 3 = ^^-[zzcoshv(L — a)sinhv(L — a) 

— vcos u[L — s)sinz<(L — a)] 


(3a) 




- A, 



c = z/vzzsinhv(L — a) sinhzz(L — a) 

— 2 l[coshv(L — a) coshzz(L — a) — l] 
c 1 = ?/[coshv(L — a) coshzz(L — a) — l] 

— 2 v /2 sinh v(L — a) sinhrz(L — a) 

c 2 = (v 2 — zz 2 )[vsinhv(L — a) coshzz(L — a) 

— zz cosh v(L — a) sinhzz(L — a)] 


v zz- , 


zzcoshv(L — a) sinhzz(L — a) 


( 2 


vcoshw(L — a) sinh v(L — a) 


(3b) 


Energy release rate failure criterion 


The deflection (5 P at the loading point may be written 

dp = dp c + dp w + Ae (4) 

where d Pc is the contribution from the cantilever according 
to ordinary beam theory 


4c = 2 


' a 3 6 a 
3EI 5GA 


(5) 


and (d Pw + d P0 ) is the contribution from the beam on elastic 
foundation 


2 1 

4 W + 4 e = 2 (w 0 + a@„) = —-(c 3 + 2 c 3 a + c 2 a 2 )/ > ( 6 ) 

El c' ’ 

The compliance, C, is given by 



2 

El 


a 6 El 1/ 2 

— +-a + — c, + 2c,a + c 7 a 

3 5 GA c 1 


)) 


(7) 


For a linear elastic body loaded by a single load, P, the crack 
propagation energy release rate, Gj , is given by 7 


P 2 dC 
2b da 


( 8 ) 


A crack starts propagating when the energy release rate 
assumes a critical value, i.e., the failure criterion is 


Q=Q, (9) 

Assuming static or quasi-static conditions and no energy 
dissipation outside the fracture region, the critical energy 
release rate is equal to the material property fracture en¬ 
ergy, Q t , i.e., 

& = 5f (io) 

From Eqs. 8-10 it follows 


Case II: X < l l i rj 2 
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P .= 


1 2 bg t 

dC 
d a 


( 11 ) 


For an infinitely long beam [(L — a) —> oo] 9 Eq. 11 leads 
to 




Jfl = /lj—— ■ 1 E 


10 G \ 3/tK 

( 12 ) 


Maximum stress failure criterion 


The maximum tensile stress in the foundation occurs at 
x = 0 and is given by 

a 0 = Kw 0 (13) 

Failure is here assumed to occur when the maximum tensile 
stress equals the perpendicular-to-grain tensile strength, 
i-e., 

o 0 =f t (14) 

From Eqs. 1, 13, and 14 it follows that 


stress, /„ and fracture energy dissipation, Gj ( . Because the 
special fracture layer has no physical perpendicular-to-grain 
dimension, strain is not defined and the perpendicular-to- 
grain stress, a, is used as a constitutive relation as a function 
of the perpendicular-to-grain deformation, d,. The fracture 
energy of the fracture layer is in general given as 


St = Jade5 f (17) 

0 

From Fig. 4, it follows that the stiffness of the linear fracture 
layer, K t , is given by 


f 2 

K f = 

2St 


(18) 


The elastic perpendicular-to-grain strain, e s , in the beam 
is given by Flook’s law, and the deformation, <5 S > and stiff¬ 
ness, K s , of the part of the DCB specimen between the two 
beam axes due to elastic strain is thus given by 


4 = he s =^rO => K s =^- (19) 

E y h 

From Fig. 3, it follows for the foundation 


P e =ft 


El c/c, 

K a + c-i/c x 


(15) 


<5(x) = 2w(x) = 4( x ) + 4(x) = 



( 20 ) 


For an infinitely long beam [(L — a) —» o°], Eq. 15 leads 
to 


h 


= bJ—hE—, ., 

'12 K{ a + Aa 


A 1 I 1 E A 

Aa = hJ -+ ,, 

1 10 G K 3 hK 

(16) 


Foundation properties 

The deformation of the foundation may be assumed to be 
composed of a contribution from a special fracture layer 
(with no physical thickness) and a contribution from the 
perpendicular-to-grain elastic strains in the beam as indi¬ 
cated in Fig. 3. 

The damage and fracture performance of wood is in 
general nonlinear, but is in the present analysis represented 
by a linear response that is equivalent in terms of peak 


or 

<5(x) = Ahv(x), 


/. 2 


K t = 

2§ f 


K= 2K t K s 


K t + K s 


E v 

= — 
h 


( 21 ) 


In previous applications of beam on elastic foundation 
theory to fracture problems, 1-5,8 the special fracture layer as 
considered in the present article has been omitted, and the 
foundation deformations have been attributed solely to the 
elastic perpendicular-to-grain strains, i.e., K t —» °° resulting 
in K = 2EJh has been assumed. If disregarding the elastic 
perpendicular-to-grain strains, i.e., K s —» °°, then 

K=2Kf=^r (22) 

yf 

From Eqs. 12 and 16 it follows immediately, that the energy 
release rate criterion and the maximum stress failure crite- 




Fig. 3. Foundation model 


Fig. 4. Constitutive relation of fracture layer 
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rion lead to the same solution if (and only if) K = f 2 IG t . This 
can also be shown to hold true for the general solutions 
given by Eqs. 11 and 15. 

The present strength analysis of mode I failure using 
beam on elastic foundation theory and introducing a special 
deformation layer with stiffness as given by Eq. 22 is a 
complete analogy to the fracture mechanics application of 
the Volkersen model to strength analysis of mode II failure 
in lap joints. 9 The analysis has been termed quasi-non-linear 
fracture mechanics because the material responses are as¬ 
sumed to be linear as in linear elastic fracture mechanics 
(LEFM). However, at the same time the tensile strength is 
assigned a finite value, not an infinite value as in LEFM, and 
a finite, nonzero size of the fracture region is considered, 
leading to a failure load that is not proportional to the 
square root of the fracture energy. 


Fracture energy 


Using the foundation modulus as given by Eq. 22, Eqs. 11 
and 15 give the same failure load and the fracture energy 
may be calculated from 


Gt = 


bEI 


Aa) 2 


P J 


P = 




Aa = -+ 


(23) 


where the c functions are given in Eq. 3, and I = bh 3 /12. 

For an infinitely long specimen \{L - a) —» oo], Eqs. 12 
and 21 lead to 


bEI 

pL 

6 h 


Gt = + Aa) 2 , Aa = h\~ + e. 


10 G 


f hE ^ 2E^ 

V Ey 


f 2 


Especially 


(24) 


Gt 




(27) 


Equations 23 and 24 do, in general, not provide explicit 
expressions for calculation of the fracture energy. However, 
a reasonable initial guess of (g t used in the right side of 
Eq. 23 or 24 usually results in convergence of (g t after just 
two to three successive calculations. Another disadvantage 
of using K =f 2 J (j, is that it involves the perpendicular-to- 
grain tensile strength of wood, which is highly volume de¬ 
pendent. A relevant size of a tension test specimen for 
direct determination of f t is not obvious. 

So-called plate joint test specimens have previously 
been proposed for deriving relevant fracture properties 
for use in LEFM models, 410 and may in principle also be 
used for determination of the perpendicular-to-grain tensile 
strength by means of a quasi-non-linear fracture mechanics 
model. 11 However, the plate joint specimen is relatively 
insensitive to variations in f in the relevant range, and 
therefore seems inappropriate. Fortunately, it appears 
that the DCB specimen is also not sensitive to variations in 
f in the relevant range for sufficiently large alh ratios. 
Figure 5 shows the failure load as predicted by Eq. 12 or 
16 and Eq. 22 using b = 20 mm, h = 50 mm, E = 8340 MPa, 
G = El 18, and = 0.23 N/mm. It may even be considered 
to use DCB specimens with large alh ratios for determina¬ 
tion of Gj t , and DCB specimens with small alh ratios for 
determination of f. 


Experimental 

Tests were conducted on DCB specimens in order to evalu¬ 
ate the difference in fracture energy determined by differ¬ 
ent theoretical models and to evaluate the influence of 
initial crack length, a , and the length of initially noncracked 
length, L — a. In addition, tests were conducted to evaluate 
the influence of elastic perpendicular-to-grain strains. 


- for /, —»t 

6 E y 

j-JW ,m e ^ 


r 2 E w 

i.e. K —» —- 


f rl\ 

i.e. K^> — 


(25) 


Gustafsson and Larsen 4 presented a solution based on 
LEFM (but not on beam on elastic foundation theory) and 
arrived at Eq. 24 with £ = 0. That solution is obtained from 
Eq. 12 if assuming K —» °°. 

If the cantilevers of the DCB specimen are considered to 
be rigidly clamped at the end of the initial crack tip, the 
compliance method leads to 


& = 


bEI 


i 2 + — -h 2 
10 G 


(26) 


ASTM D3433-93 5 gives the following formula for determi¬ 
nation of the fracture energy 



/,[ MPa] 

Fig. 5. Influence of f on the failure load 
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Methods and materials 

All specimens were made of glulam of Japanese cedar 
('Cryptomeria japonica). The laminae, thickness 30mm, 
were all of the same grade and without finger joints. The 
MOE in the grain direction was determined by measuring 
the longitudinal vibration frequency of the glulam beams, 
from which the specimens were cut. All tests were displace¬ 
ment controlled. All DCB specimens were tested so that 
failure occurred after 2-3 min, the time to failure for the 
tension specimens was about 30 seconds. 

Three test series were conducted on DCB specimens, as 
shown in Fig. 1. The initial crack was made by drilling a 
3-mm hole at the location of the initial crack tip and then 
cutting a 2-mm-wide slit from the end of the specimen to the 
hole. 

Series 1 was conducted to evaluate the influence of the 
initial crack length, a. All specimens had the dimensions of 
b = 20mm, h = 50mm, L = 600mm. Ratios alh = 1, 2, 4, 
and 6 were tested. Five specimens were tested for each 
condition. Moisture content (MC) was 11%, the density was 
370 kg/m 3 at 11% MC, and the MOE was 8340 MPa. The 
specimens were not specially selected to avoid defects in the 
vicinity of the crack tip. 

Series 2 was conducted to evaluate the influence of the 
length of the initially uncracked part, L — a. All specimens 
had the dimensions of b = 20 mm, h = 50 mm, a = 100 mm. 
The number of specimens were 6, 12, 10, and 12 for the 
tested ratios (L — a)lh = 2, 4, 6, and 12, respectively. 
The MC was 14%, density was 384kg/m 3 at 14% MC, and 
the MOE was 8959 MPa. 

Series 3 was conducted to evaluate the influence of beam 
depth, h. Specimens with the dimensions of b = 20mm, 
a = 2 h, and (L — a) = Ah were tested for h = 20mm, 50mm, 
and 100 mm. Six specimens were tested for each condition. 
The MC was 13%, density was 358 kg/m 3 at 13% MC, and 
the MOE was 7914 MPa. All specimens for this test series 
were cut from the same glulam beam, and in such a way that 
the initial cracks were all located in the same lamina. The 
specimens were carefully selected to avoid defects in the 
vicinity of the crack tip. 

Tension tests were conducted on specimens as shown in 
Fig. 6. The specimens were cut from the same beam as the 
DCB specimens of test series 2, and in such a way that 
failure in the two specimen types occurred in the same 
lamina. Six specimens were tested. 


Fig. 6. Test specimen for determination of perpendicular-to-grain 
tensile strength 


Results and discussion 

The tension tests gave a value of the perpendicular-to-grain 
tensile strength,/,, of 3.9 ± 0.9MPa (mean value ± standard 
deviation). Table 1 shows the mean values of the fracture 
energy as determined by the use of Eqs. 24,26, and 27 given 
for test series 1. ft may be worthy of notice that the more 
elaborate a model (i.e., the more contributions are taken 
into account for calculation of the compliance), the higher 
values of the fracture energy are obtained. This fact may be 
a general point of concern. If for instance Eq. 24 is used in 
a testing standard for determination of the fracture energy, 
and this value of the fracture energy is used for calculation 
of the load-carrying capacity using a less elaborate model 
(say Eq. 27), an unsafe solution is obtained. Using K = f 2 J(g f 
will always produce safe values of the fracture energy as 
compared with K as given by Eq. 21, and in most practical 
cases will also produce safe values as compared with 
K = 2 E y lh. 

Figure 7 shows the failure loads of test series 2. The mean 
values of the failure loads for (L — a)lh = 4, 6, and 12 are 
not significantly different, while (L — a)/h = 2 gives a 
significantly lower mean failure load. The theoretical failure 
load as predicted by Eqs. 15 and 22 using (g f = 0.23 N/mm 
(which is the mean value determined from test series 1 for al 
h = 2,4, and 6) is also shown for f = 3.9 MPa as obtained by 
the tension tests and for a value f = 2.3 MPa chosen in order 


Table 1. Fracture energy determined for test series 1 


Equation 

alh 




1 

2 

4 

6 

24 

0.26 

0.27 

0.27 

0.24 

24 (f t -> °o) 

0.25 

0.26 

0.26 

0.24 

24 ( E y oo) 

0.20 

0.22 

0.23 

0.22 

24 (s = 0) 

0.16 

0.19 

0.21 

0.20 

26 

0.08 

0.10 

0.13 

0.14 

27 

0.04 

0.07 

0.12 

0.14 


Units: N/mm 



Fig. 7. Failure loads of test series 2 
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Table 2. Fracture energy of test series 3 using Eq. 24 


Specimen 

h (mm) 

a (mm) 

Gj[ (N/mm) 

E y -> °° 

8 

T 

1 

20 

40 

0.16 

0.18 

2 

50 

100 

0.15 

0.18 

3 

100 

200 

0.18 

0.23 


Table 3. Calculated and tested failure loads of test series 3 

Equation 

K 

Pa (N) 

T c3 (N) 

PJP, 3 

12 


237 

558 

0.42 

16 

fJQ i 

237 

558 

0.42 

12 

2E y lli 

244 

546 

0.45 

16 

2 E y lh 

402 

2010 

0.20 

Test 


248 

614 

0.40 


P cl , failure load of specimen 1 (h = 20 mm ); p*, failure load of speci¬ 
men 3 (h = 100 mm) 


to fit the test results well. The model is seen to be able 
to predict the influence of the initially uncracked length, 
L — a, in good agreement with the test results. 

Table 2 shows the fracture energy determined for series 
3 using Eq. 24 with E y —» °° (i.e., K = fj(g^ and/, —> °° (i.e., 
K = 2E y /h). The fracture energy of test series 3 is 
significantly lower than the values determined for test series 
1 and 2. This is due to the fact that all specimens of test 
series 3 were taken from the same beam and carefully se¬ 
lected so as to avoid any defects, which have a reinforcing 
effect, in the vicinity of the crack tip. 

Table 3 shows the failure loads, P cl and P«, of test speci¬ 
mens 1 and 3 (see Table 2), respectively, as obtained by 
testing and as predicted by Eqs. 12 and 16 using K = f 2 l<g f 
(i.e., E y —> °°) and K = 2 E y lh (i.e.,/, —> °°). Specimen 2 ( h = 
50 mm, a = 100 mm) was used for determination of the 
fracture energy as given in Table 2. Because a constant 
value alh = 2 was used for all specimens, Eq. 12 with K — 
2E y lh leads to //,///, = (/;,/ h 3 ) 112 and Eq. 16 with K = 2 E y l 
h leads to P cl IP c3 = hj h 3 , which are independent of/ and Gj,. 

Table 3 also shows that the maximum stress failure crite¬ 
rion (Eq. 16) is in severe disagreement with test results if 
using K = 2E.jh, while K = f\t(g t leads to agreement. How¬ 
ever, the two different foundation moduli both result in 
good agreement between theory and tests if used with the 
compliance approach (Eq. 12). 


Conclusions 

A model for analysis of the double cantilever beam speci¬ 
men for determination of the fracture energy of wood was 
presented. The model is based on linear elastic Timoshenko 
beam on elastic foundation theory, which includes shear 
deformations and takes into account the finite length of the 
initially uncracked part of the DCB specimen. The analysis, 
which may be characterized as a quasi-non-lincar fracture 
mechanics analysis, includes a fracture mechanics approach 


based on the compliance method and a conventional stress 
approach. It was shown that a particular choice of the prop¬ 
erties of the elastic foundation, which ensures that the 
perpendicular-to-grain tensile strength and fracture energy 
properties of the wood are correctly represented, makes the 
two approaches lead to the same solution. Previous models 
using beam on elastic foundation theory seem all to have 
related the stiffness of the foundation solely to the elastic 
perpendicular-to-grain strains, which leads to the fact that 
the compliance method and stress approaches, in general, 
lead to significantly different results. 

Tests were conducted to evaluate the influence of the 
length of the initial crack and of the length of the initially 
uncracked (beam on elastic foundation) part of the DCB 
specimen. Theoretically predicted results were found to be 
in very good agreement with test results. 

Tests were also conducted on DCB specimens with vari¬ 
ous depths in order to evaluate the appropriate choice of 
foundation properties. The tests revealed that relating the 
foundation stiffness to the elastic perpendicular-to-grain 
strains leads to good agreement between theory and tests if 
used with the compliance method approach, but to very 
poor results if used with the stress approach. The new esti¬ 
mation of the foundation stiffness as presented in this 
article unifies the compliance method and stress ap¬ 
proaches and leads to good agreement with test results. 
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